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Abstract 



We report on a search for i?-parity violating super symmetry in pp collisions 
at = 1.8 TeV using the D0 detector at Fermilab. Events with at least 
two electrons and four or more jets were studied. We observe 2 events in 
99 ± 4.4 pb _1 of data, consistent with the expected background of 1.8 ± 0.4 
events. This result is interpreted within the framework of minimal low-energy 
supergravity supersymmetry models. Squarks with mass below 243 GeV/c 2 
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and gluinos with mass below 227 GeV/c 2 are excluded at the 95% confidence 
level (C. L.) for Aq = 0, \i < 0, tan (3 = 2 and a finite value for any one of 
the six i?-parity violating couplings \\jk (j =1, 2 and k =1, 2, 3). 

The standard model (SM) has survived many precision tests. However, it is thought 
incomplete, and supersymmetry (SUSY) |l| is considered an attractive extension to the 
SM because it protects the Higgs mass from large radiative corrections and can provide a 
dynamical means for breaking electroweak symmetry. SUSY predicts for each particle in 
the SM, a partner with spin differing by half a unit. In its general form, the theory contains 
over one hundred free parameters. For our comparison with data, we have therefore chosen 
the more tractable framework provided by minimal low-energy supergravity (mSUGRA) , 
which has only five free parameters: a common mass for scalars (mo), a common mass for all 
gauginos (rrii/z), and a common trilinear coupling constant (Aq), all specified at the grand 
unification scale. The other two parameters are the ratio of the vacuum expectation values 
of the two Higgs doublets (tan j3), and the sign of the Higgsino mass parameter /x. The 
masses and couplings at the weak scale are obtained from these five parameters by solving 
a set of renormalization group equations. 

Most of the searches for supersymmetric particles reported thus far have assumed the 
conservation of a multiplicative quantum number called i?-parity ||. i?-parity is defined as 
R = (— i) 3B + L + 2S ; where B, L and S are the baryon, lepton and spin quantum numbers, 
respectively. R is +1 for SM particles, and —1 for their SUSY partners. In SUSY, f?-parity 
violation can occur quite naturally through the following Yukawa coupling terms in the 
superpotential: 

KjkLiLjEk + \j k LiQjDk + \j k UiDjDk 

where L and Q are the SU(2)-doublet lepton and quark superfields; E, U, and D are the 
singlet lepton, up and down type quark superfields, respectively; and i, j and k are the 
generation indices. The Yukawa couplings are antisymmetric in the same superfield indices. 
Thus, there can be up to 45 new Yukawa terms. We have therefore made the following 
simplifying assumptions for our analysis: 

• Among the 45 i?-parity violating couplings, only one dominates. This is motivated by 
the fact that the new couplings are similar to the SM Yukawa couplings, for which the 
top quark Yukawa term dominates. Moreover, bounds on products of two couplings 
are generally stringent, because the presence of more than one coupling can induce 
rare processes such as flavor changing neutral currents at the tree level Qj. 

• The strength of the .R-parity violating coupling under consideration is > 10~ 3 , so that 
the lightest supersymmetric particle (LSP) decays close to the interaction vertex. This 
is consistent with the existing upper bounds on the strength of the couplings from low 
energy experiments ||. 

• The strength of the finite imparity violating coupling term is significantly smaller 
than the gauge couplings. Thus, supersymmetric particles are produced in pairs, and 
i?-parity violation manifests itself only in the decay of the LSP. 
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Of the three kinds of Yukawa coupling terms, the 5-violating A" are difficult to study 
at the Fermilab Tevatron as they lead to events with multiple jets that would be over- 
whelmed by large backgrounds from QCD production of jets. However, the L-violating A 
and A' couplings give rise to multilepton and multijet final states [|], which provide excellent 
signatures. 

This Letter reports on an analysis of the dielectron and four jets channel, interpreted 
in the mSUGRA framework, with i?-parity violating decays of the LSP. In mSUGRA, the 
lightest neutralino is almost always the LSP except in a small region of the (mo, 7711/2) plane 
where the sneutrino is the LSP (indicated in Fig. [I]). But the mass of the sneutrino in that 
region is below 39 GeV/c 2 and, hence, excluded (> 43.1 GeV/c 2 at 95% C. L.) [0 by the 
known invisible decay width of the Z boson, assuming that there are three degenerate left 
handed sneutrino species. We assume that all the .R-parity violating couplings are small 
except for one of the six A' ljfc (j =1, 2 and k =1, 2, 3), so that each LSP decays into one 
electron and two quarks which gives rise to final states with two or more electrons and four 
or more jets that we consider in our analysis. 

The D0 detector || has three major subsystems: a central tracker, an uranium liquid 
argon sampling calorimeter, and a muon spectrometer. Electrons are identified as narrow 
energy clusters that deposit more than 90% of their energy in the electromagnetic sections 
of the calorimeter. Jets are reconstructed using a cone algorithm || with radius 0.5 in 
pseudorapidity — azimuthal angle (77, (f>) space. The data used for this analysis were collected 
during the 1994-1995 Fermilab Tevatron run at a center-of-mass energy of 1.8 TeV, and 
correspond to an integrated luminosity of 99 ± 4.4 pb _1 []T0 . 

Our initial sample of 163,140 events was collected with triggers requiring at least five 
calorimeter energy clusters, and Ht > 115 GeV, where Ht is the scalar sum of the transverse 
energies (Et) of all calorimeter clusters. In the offline analysis, we required at least two 
electrons, one with Et > 15 GeV and the second with Et > 10 GeV, and at least four 
jets with Et > 15 GeV. Electrons had to be either within \r/\ < 1.1 (central calorimeter) or 
1.5 < \r/\ < 2.5 (forward calorimeters), to be isolated from other energy deposits, and to have 
shower shape and tracking information consistent with that expected for electrons ]TT], 12] . 



Jets had to be within \r]\ < 2.5. The requirement on electrons reduced the original sample 
to just 38 events, and the subsequent requirements on jets reduced it further to 6 events. To 
suppress backgrounds from electron decays of Z bosons, we rejected events whose dielectron 
invariant mass was in the range of 76-106 GeV/c 2 . To ensure high trigger efficiency, events 
were further required to have Ht > 150 GeV. The cut on Z mass reduced our data sample 
to 2 events, but the Ht requirement had no further impact. 

The major inherent SM backgrounds are from Drell-Yan production (DY), from the 
decay of tt to electrons, and from the decay of Z bosons to r pairs that subsequently decay 
to electrons. Events arising from the misidentification of jets as electrons comprise the major 
source of instrumental background for this analysis. The huge reduction in our data sample 
from the requirement of having two isolated electrons reflects the fact that most of the events 
passing the trigger are due to QCD multijet production, and have no true isolated electrons. 

A geant |13| based simulation of the D0 detector was used to estimate efficiencies 
of the kinematic cuts for non-instrumental backgrounds. Measured electron identification 
efficiencies were then folded in to calculate the net detection efficiency. Using Z(— > ee) + jets 
data, we estimated single-electron identification efficiencies to be 0.68 ± 0.07 in the central 
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calorimeter, and 0.60 ± 0.07 in the forward calorimeters, isajet [|14| was used to generate 
DY events, with cross section increased by a factor of 1.7 to obtain agreement with the Z 
+ multijet data in the mass region of the Z boson, yielding an expected 0.37 ± 0.14 (stat.) 



±0.14 (syst.) events. Top quark events were generated using the herwig [15] program. 
The measured cross section for ti production (5.9 ± 1.7 pb) [|16j was used to estimate this 
contribution to background to be 0.07 ± 0.02 ± 0.02 events. The production cross section 
of the Z boson multiplied by its leptonic branching fraction of (221 ±11) pb [[TO] was used 
to estimate the background due to Z(— > rr — > ee) to be 0.07 ± 0.01 ± 0.02 events. The 
instrumental background was estimated from data in two steps. First, from multijet data, 
we estimated the probability for misidentifying a jet as an isolated electron. This was (4.6 ± 
0.4) x 10~ 4 in the central and (1.4±0.2) x 10~ 3 in the forward calorimeters. Within statistical 
accuracy, these probabilities were found to be independent of Et- We then selected a multijet 
data sample passing the same kinematic requirements as our data sample, but requiring two 
additional jets instead of two electrons. The number of background events was estimated 
to be 1.27 ± 0.24 (with negligible statistical uncertainty) by applying the probability for jet 
misidentification to these multijet data. The statistical components of uncertainty include 
fluctuations due to the finite sample size of simulated events and uncertainties in electron 
identification efficiencies. The systematic components of the uncertainty include those due 
to jet energy scale and values of cross sections. Our two observed events are consistent with 
the expected background, both in the number of expected events 1.8 ± 0.2 ± 0.3, and in 
their kinematic characteristics. In what follows, we interpret this null result in terms of an 
excluded region in mSUGRA parameter space. 

Using isajet, we generated signal events at 125 points in the (m^m^) plane, with 
A = 0, < and tan (3 = 2. i?-parity violating decays of the LSP are not available in 
isajet. The desired decay modes and branching fractions for the LSP were therefore added 
separately. The branching fraction of the LSP into a charged lepton or neutrino depends 
on the gauge composition of the LSP, which in turn depends on the mSUGRA parameters. 
This was incorporated into isajet using the calculation of Ref. [|17]]. Once we specify a decay 
mode, isajet does a 3 body phase-space decay, but does not implement the appropriate 
matrix element into the differential distribution. The efficiency multiplied by the branching 
fraction for each signal sample was determined using a method similar to that used for the 
estimation of the SM background. The expected event yields in the (m ,mi/ 2 ) parameter 
space, corresponding to our integrated luminosity of 99 pb -1 , are given in Table |. 

For each point in the (roo,™^) plane, we obtained a 95% C. L. upper limit on the 
cross section for signal. This was done using a Bayesian technique, with a flat prior for 
the signal cross section, and Gaussian priors for the luminosity, efficiency, and expected 
background. The excluded region in the (mo, mi/2) plane was then obtained by comparing 
the limits on the measured cross section with the leading-order SUSY prediction given by 
isajet. This is shown in Fig. The slanted hatched area in Fig. 1 indicates the region 
in which the model does not produce radiative electroweak symmetry breaking. In the low 
m region (m < 150 GeV/c 2 ), the dominant SUSY process that contributes to the signal 
is pair production of squarks. Hence, in this region, the exclusion contour follows a squark 
mass contour (m^ = 273 GeV/c 2 ). The dip in the contour for mo = 60 — 80 GeV/c 2 can 
be attributed to the fact that the two electrons can originate either from the decay of LSPs 
or from other SUSY particles. In about 60% of the cases, both LSPs decay into electrons. 
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TABLE I. Efficiency (e) multiplied by the branching fraction (B) and the expected event yield 
(N), for several points in the (m-cb m i/2) parameter space. The uncertainties are the sum in 
quadrature of the statistical and systematic uncertainties (the statistical uncertainty dominates). 



m ( GeV/c 2 ) 


m 1/2 ( GeV/c 2 ) 


eB{%) 


(AT) 





120 


1.59 ±0.23 


3.5 ±0.5 


50 


110 


1.49 ±0.22 


2.8 ±0.4 


120 


110 


1.86 ±0.25 


3.3 ±0.4 


190 


100 


1.56 ±0.22 


3.4 ±0.4 


280 


90 


0.95 ±0.15 


2.9 ±0.4 


320 


90 


0.71 ±0.13 


2.2 ±0.4 



However, electrons arising from the decay of LSPs may not always pass the cut. In 
such cases additional electrons arising from the decay of the second lightest neutralino (x°) 
can make the event pass the Et criterion. But for ttiq = 60 — 80 GeV/c 2 , sneutrinos 
become lighter than the an d the decay of X2 t° Xi an d neutrinos (x° — * vv\ v ~* Xi p ) 
becomes dominant. This reduces the overall branching fraction to dielectrons, resulting in 
the observed dip. 

As mo increases, the sneutrino becomes heavier than x®, an d consequently the branching 
fraction of x\ to neutrinos decreases, leading to an increase in the rate for the competing 
selectron channel, thereby enhancing the branching into the dielectron mode. (That is, when 
the X2 decay proceeds through a virtual sneutrino, the decay through a virtual selectron 
becomes competitive.) The exclusion contour therefore moves up and again follows the 
273 GeV/c 2 squark mass curve until the intermediate m region (150 GeV/c 2 < m < 
280 GeV/c 2 ), where processes such as the production of gluinos, Xi-i an d X2, start becoming 
important. The masses of these particles, as well as their production cross sections, do not 
change much with the increase of mo- As a result, the exclusion contour in this region 
becomes less dependent on m . 

Finally, in the asymptotic region (m > 280 GeV/c 2 ), production of squarks becomes 
insignificant, and the contour of exclusion becomes totally independent on mo- In Fig. |], 
we have overlaid contours of fixed gluino mass and the average of the masses of the first 
two generations of squarks. Squarks with mass below 243 GeV/c 2 and gluinos below 227 
GeV/c 2 are excluded for A = 0, \i < 0, tan (3 = 2, and a finite value (> 10~ 3 ) for any 
one of the six A' ljfc (j =1, 2 and k =1, 2, 3) couplings. For equal mass squarks and gluinos, 
the corresponding limit is 277 GeV/c 2 . 

We note that our results are essentially independent of the choice of A , as it affects only 
third generation sparticle masses. For [i > and higher values of tan /?, the sensitivity of our 
search is expected to fall for two reasons: 1) the photino component of the LSP decreases, 
resulting in the decrease of the branching fraction of the LSP into electrons; 2) the charginos 
and neutralinos become light, resulting in events with softer electrons and jets that fail the 
kinematic requirements. We have estimated the sensitivity of our search for larger values 
of tan P, by extrapolating our tan (3 = 2 results using smeared parton level isajet W2 



(without full detector simulation). Figure 2 shows the region excluded at 95% C. L. in the 
(m , rrii/2) plane for A = 0, /i < 0, tan (3 = 6, and a value > 10~ 3 for any one of the six 
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i i i r 




50 100 150 200 250 300 350 400 

m (GeV/c 2 ) 

FIG. 1. Exclusion contour in the (tjo> ^1/2) plane for Aq = 0, fi < 0, tan (3=2 and a finite 
X' 1 j k (j =1, 2 and k =1, 2, 3) coupling. The region below the bold line is excluded at the 95% C.L. 
The slanted hatched region is excluded for theoretical reasons. In the horizontally hatched region, 
the sneutrino is the LSP, but is excluded by searches at LEP (see the text). 



K.jk = 1' ^ and k =1, 2, 3) couplings. For higher values of tan (3, the sensitivity of this 
search deteriorates rapidly and requires a different analysis. 

In conclusion, we have searched for events containing at least two electrons and four or 
more jets. Such events would be characteristic of processes involving the pair production of 
SUSY particles with the decay of the LSP through a .R-parity violating coupling. Finding no 
excess of events beyond the prediction of the standard model, we interpret this result within 
the mSUGRA framework as an excluded region in the (m , mi/2) plane for fixed values of A 
and sign of // and for several values of tan f3. This is the first result reported from Tevatron 
on search for .R-parity violating SUSY involving several A' couplings in the mSUGRA frame- 
work. The Tevatron will continue to provide unique opportunities for searching for f?-parity 
violating SUSY in a larger range of parameter space with the improved data anticipated 



from the next run 12 
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50 100 150 200 250 300 350 400 

m (GeV/c 2 ) 

FIG. 2. Exclusion contour in the (mo, 011/2) plane for ^0 = 0,^ < 0, tan (3 = 6, and a finite 
Kjk ( J =1> 2 an d fc =1, 2, 3) coupling. 
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